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NACA RM L 54 A 20 CONFIDENTIAL 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

RESEARCH MEMORANDUM 

AERODYNAMIC CHARACTERISTICS AT TRANSONIC AND SUPERSONIC 
SPEEDS OF A ROCKET-PROPELLED AIRPLANE CONFIGURATION 
HAVING A 52 . 5° DELTA WING AND A LOW, 

SWEPT HORIZONTAL TAIL 
By Alan B. Kehlet 

SUMMARY 


A flight investigation over a Mach number range from 0-79 to I.83 
has been conducted in order to determine the aerodynamic characteristics 
at low lift of a rocket model of an airplane configuration having a 
52 . 5° delta wing of aspect ratio 3.08 with NACA 65AOO3 airfoil sections 
in the streamwise direction and a low, swept horizontal tail. The 
lift-curve slopes and static longitudinal stability were nonlinear with 
lift coefficient over most of the Mach number range and increased with 
increasing lift coefficient over the lift range covered. The minimum 
drag coefficient increased from 0.016 to 0.035 through the transonic 
speed range. 

The model damping characteristics were irregular and altered for 
the two tail settings over the Mach number range covered with the higher- 
lift tail setting having the greater damping. 

The measured periods of the lateral oscillations were of the same 
order of magnitude as the longitudinal periods, but apparently were not 
affected by lift coefficient. The model exhibited stable static direc- 
tional characteristics throughout the Mach number range tested. 


INTRODUCTION 


As part of a general research program investigating longitudinal 
stability of wings having various plan forms and thickness ratios, a 
rocket -propelled model of an airplane configuration having a 52. 5° delta 
wing of aspect ratio 3-08 has been flown. The basic fuselage-empennage 
configuration (ref. l) had swept horizontal and vertical tails with the 
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all-movable horizontal tail mounted in a low position. During the 
flight , the horizontal tail was deflected in a square-wave program 
between 0.10° and -3.2°. 

The model was flown at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 
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SYMBOLS 

normal-force coefficient, 

8 q 

chord-force coefficient, — - 

g q 

lift coefficient, C]y cos a - Cg sin a 

drag coefficient, Cg cos a + sin a 

lift coefficient at minimum drag coefficient 

pitching-moment coefficient about center of gravity 

normal acceleration determined from accelerometer, ft/sec^ 

longitudinal acceleration determined from accelerometer, 
ft/sec^ 

acceleration of gravity, ft/sec^ 
dynamic pressure, 0 . 70 pM^ 
free-stream static pressure, lb/sq ft 
Mach number 

Reynolds number, based on wing mean aerodynamic chord 
wing area (including area enclosed by fuselage), sq ft 
wing mean aerodynamic chord (M.A.C.), ft 
wing span, ft 
weight, lb 
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cross-sectional area, sq in. 

equivalent radius of cross-sectional area, in. 

longitudinal distance from station 0 , in. 

length of model, in. 

period of oscillation, sec 

time to damp to one-half amplitude, sec 

angle of attack, deg 

control panel deflection (measured free stream normal to 
wing chord plane), deg 

angle of pitch, radians 

angle of sideslip, deg 


dC m 




Crn d " 


<*C m 



c 


n |3* 


effective rate of change of yawing-moment coefficient with 
sideslip angle per degree (derived as in ref. 7 ) 


dc D = ac D 
ac L 2 d (c L - c^,) 2 

Subscripts: 
w wing 

f fuselage 
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1 da 

57-3 dt 


q. = 


d9 

dt 


The symbols a, 
the derivative of the 


example, 



a, q, 5, and p* used as subscripts indicate 
quantity with respect to the subscript; for 


MODEL AND INSTRUMENTATION 
Model 


A three-view drawing of the model is shown in figure 1 . The non- 
dimensional equivalent body and area distribution, presented for 
possible drag correlation purposes, are shown in figure 2 . Photographs 
of the model are shown in figure 3 * 

The empennage has a vertical fin of wood and aluminum with the 
quarter-chord line swept 60 ° and NACA 65 AOO 3 airfoil sections in the 
streamwise direction; a horizontal tail of duralumin with the quarter- 
chord line swept 45° with 20° negative dihedral and NACA 65 AOO 6 airfoil 
sections in the streamwise direction perpendicular to the chord plane. 
The fuselage is described in reference 2. 

The steel delta wing of aspect ratio 3.08 had a leading-edge sweep 
of 52 . 5 ° NACA 65 AOO 3 airfoil sections in the streamwise direction. 

Each panel of the horizontal tail was deflected in an approximate 
square-wave program by a separate servo-control fed by a common pressure 
system and regulated by an electric motor-driven selector valve. For 
the present investigation, the stop positions were 0 . 1 ° and - 3 . 2 ° 
measured parallel to the model center line and normal to the wing chord 
plane. 

The model weighed 118.75 pounds and had a moment of inertia in 
pitch and yaw of 8.21 and 8-31 slug- ft 2 , respectively. The center of 
gravity was located at 0.2 6 of the wing mean aerodynamic chord. 
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Instrumentation 

The model was equipped with an NACA telemetering system which 
transmitted continuous measurements of normal acceleration at the 
center of gravity, normal acceleration at a reference nose station, 
angle of attack, longitudinal acceleration, transverse acceleration, 
control position, total pressure, and reference static pressure. 

Flight-path information was obtained from tracking radar and 
atmospheric conditions at altitude from a radiosonde released immediately 
after the flight. Motion-picture cameras were used to photograph the 
launching and first portion of the flight. 


TESTS AND DATA REDUCTION 
Preflight Tests 

Prior to instrumentation, the model was placed in a profile machine 
and measurements were made in the streamwise direction of the airfoil 
sections, perpendicular to the chord planes of the vertical and hori- 
zontal tails and of the wing. These measurements were used to check 
construction tolerances. 

Prior to flight testing and with the instruments installed, the 
model was suspended by shock chords and vibrated by an electromagnetic 
shaker. The following model natural frequencies and modes of vibration 
were determined: 



First 

bending, 

cps 

Second 

bending, 

cps 

Torsion, 

cps 

Horizontal tail . . . 

100 



Vertical tail .... 

57 



Wing 

116 

1 64 

370 


First bending 


Nose 



Wing Modes 
Second bending 



Torsion 
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Flight Tests 

The model was launched at an angle of approximately 6o° from the 
horizontal by means of a mobile launcher as shown in figure 3 . Two 
6-inch-diameter solid-fuel ABL Deacon rocket motors boosted the model 
to maximum velocity. All measurements used were taken during the 
decelerating portion of the flight. 


Data Reduction 

The response of the model to deflections of an all-movable hori- 
zontal tail in an approximate square- wave program was analyzed by the 
method of reference 3- The indicated angles of attack were corrected 
to angles of attack at the model center of gravity by the method of 
reference 4. The two-accelerometer method for obtaining instantaneous 
total pitching-moment coefficients was used as described in reference 2. 


ACCURACY 


The absolute accuracy of the measured quantities is impossible to 
establish because the instrument calibrations can not be checked during 
or after the flight. Most of the probable instrumentation errors occur 
as errors in absolute magnitude. Incremental values or slopes should, 
in general, be more accurate than the absolute values. An indication 
of the systematic instrument errors possible is given by the following 
table, based on an accuracy of ±1 percent of the full instrument range: 


M 

C N 

Cc 

1-7 

1.0 

.8 

± 0.00490 

±.02251 

±.o4l8o 

±0.00122 

±.00563 

±.01045 


The CW Doppler radar unit is believed to be accurate to better 
than 1 precent for nonmaneuvering models. The Mach number at peak 
velocity should, therefore, be accurate to 1 percent or better. Since 
the Mach number subsequent to peak velocity was determined from 
telemetric data, it probably becomes less accurate as the Mach number 
decreases . 

Further errors in the aerodynamic coefficients may arise from 
possible dynamic-pressure inaccuracies which are approximately twice 
as great as errors in Mach number. 
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An indication of random errors encountered may tie noted from the 
scatter of data points shown in the figures. Errors in angle of attack 
and control panel deflection are independent of dynamic pressure and 
are not likely to vary with Mach number. The horizontal- tail deflec- 
tions are estimated to be accurate within ±0.10° and the increments in 
angle of attack within ±0.20°. 


RESULTS AND DISCUSSION 


Dynamic pressure and Reynolds number obtained during the flight 
are shown in figures b and 5 > respectively. A typical time history at 
the higher Mach numbers of some of the quantities measured is shown in 
figure 6. Hereafter, in order to distinguish the two tail settings, 
the deflection of 0.1° shall be referred to as the low- lift tail setting 
and the deflection of -3.2° as the higher-lift tail setting. 


Longitudinal Trim 

The variations of the trim lift coefficient and trim angle of attack 
at the two control settings as functions of Mach number are shown in fig- 
ure 7* At transonic speeds and with increasing Mach number, the model 
exhibited a trim change of approximately 1° nose-up in the low lift 
range; whereas at the higher trimmed lift condition, a nose-down trim 
change of approximately 2° occurred. 


Lift 

The variation of the lift coefficient with angle of attack over the 
Mach number range is shown in figure 8. The lift-curve slopes repre- 
sented by the faired lines in figure 8 are presented as functions of 
Mach number in figure 9 for' both tail settings. Included in figure 9 
for comparative purposes are the lift-curve-slope data for a wind-tunnel 
model with a similar wing- fuselage combination having no horizontal tail 
(ref. 5 )• I n view of the differences in models, the agreement between 
the rocket-model data and the wind-tunnel data is considered good. It 
is indicated in figure 9 that the lift coefficient does not vary line- 
arly with the angle of attack, particularly at the lower Mach numbers; 
however, reference to the trim lift coefficients (fig. 7), which are 
indicative of the lift ranges covered by the two tail settings, shows 
that the linearity of the lift-curve slopes present at the higher Mach 
numbers may be due to the low lift range of the higher-lift tail setting 
and not entirely to a Mach number effect . The increase in lift-curve 
slope with increasing lift coefficient, where nonlinearity occurs, is 
believed to be due in part to a decrease in the downwash field over the 
low horizontal tail with increasing angle of attack of the model and to 
a nonlinearity in the lift-curve slope of the long nose section of the 
fuselage . 
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Drag 

The variation of drag coefficient with lift coefficient corresponding 
to the lift ranges of figure 8 is shown in figure 10. The maximum lift- 
drag ratios that could be measured from figure 10 and lift coefficient 
at which (L/D) ma x occurs are shown in figure 11. Because of the 
reduced amplitude of the oscillations , the model never reached its maxi- 
mum lift-drag ratios at Mach numbers greater than 1.12; values of 
(L/D) m ax Cl for (k/D) max obtained at higher Mach numbers were 

extrapolated from the data for the drag due to lift, minimum drag coef- 
ficient, and lift coefficient at minimum drag. At high subsonic speeds, 
the maximum value of the maximum lift-drag ratios obtained was about 9, 
decreasing to about 6 at low supersonic speeds and to about 4^ at the 

test s limit; the values of lift coefficient corresponding to these 
values are about 0.25, 0.40, and 0.30, respectively. 

The minimum drag coefficient and the lift coefficient at minimum 
drag coefficient obtained from figure 10 are presented as a function 
of Mach number in figure 12. The flagged symbols on the minimum-drag- 
coefficient curve indicate extrapolated points from the higher- lift 
tail setting. It should be noted that, as indicated by the basic data 
of figure 10, the values of the minimum drag coefficient correspond 
closely to the values of the zero lift-drag coefficient, the drag polars 
at the low- lift tail setting are flat around zero lift coefficient, and 
the values of the lift coefficient at minimum drag can be determined, 
with any degree of accuracy, only below a Mach number of 1.3. At all 
Mach numbers for the higher-lift tail setting, no lift coefficient at 
minimum drag can be obtained. 

By using the values of the lift coefficient at minimum drag coeffi- 
cient (fig. 12) for the low-lift tail setting and by interpolating 
values for the higher-lift tail setting, the effect of lift on drag as 
a function of Mach number was determined and is presented in figure 13. 

As indicated by the data of figure 13 , the model exhibited poor leading- 
edge suction over the Mach number range where the effect of lift on drag 
could be determined. As stated before, prior to instrumentation, measure- 
ments of the airfoil sections were made. The wing airfoil sections, as 
determined from these measurements, revealed a leading-edge radius some- 
what smaller than the true NACA 65AOO3 airfoil section. The resulting 
rather sharp leading edge is believed to be part of the cause of the 
lack of leading-edge suction obtained. 


Longitudinal Static Stability 

The measured periods of the longitudinal oscillations for both 
tail settings as a function of Mach number and the longitudinal-static- 
stability parameter Cfl^ determined from these periods are presented in 
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figures 14 and 15, respectively. As with the lift-curve slopes (fig. 9), 
the nonlinearity present in the data decreased with increasing Mach num- 
ber. Where nonlinearity of the data occurred, the higher-lift tail 
setting had the greater value of static stability. Again, however, the 
linearity at the higher Mach numbers may be due to the decreasing trim 
lift coefficient rather than to the increasing Mach number. 

The variation of the static-stability parameter is reflected in 
the aerodynamic -center location for the two tail settings (fig. 1 6). 

The data at both tail settings exhibited a rearward movement of the 
aerodynamic center up to a Mach number of about 1.5 with the higher- 
lift tail setting having the greater stability. Above M = 1.5, a 
forward movement occurs . 

The variation of the total pitching-moment coefficient, obtained 
from the two- accelerometer method (ref. 2), with lift coefficient is 
shown in figure 17- Although some scatter is present, particularly at 
subsonic speeds, the data agree, in general, with the slopes indicated by 
the period method and the trim data of figure 7 which are also shown 
in figure 17. 

A measure of the horizontal-tail effectiveness in producing moment 
and ability to produce lift, as obtained from an average of Cn^ for 

the two tail settings (fig. 15), the trim angle of attack (fig. 7), and 
an estimated tail length, is presented in figure 18. Both parameters 
exhibit the same general shape as the lift-curve slopes (fig. 9), that 
is, increasing with Mach number to M = 1.0 and then decreasing with 
increasing Mach number. At the test limit (M I.83), a decrease of 
approximately 50 percent of that at M = 1.0 is noted in both param- 
eters . Included in figure 18 for comparative purposes are Cj^ and 
Cmg, obtained from from the exposed duralumin wing of reference 2 

(extrapolated Ci^ at the higher Mach numbers) and corrected for 
dihedral effects. 

Although the values obtained from reference 2 are approximately 
20 percent higher throughout the Mach number range covered, it should 
be pointed out that the duralumin wing Ci^ used was not corrected 
for the small gap that exists between the fuselage and horizontal-tail 
panel. The effects of a small gap in a nonviscous flow (ref. 6) can 
be large, reducing the lift-curve slope by as much as kO percent. By 
using an average gap of O.O76 inch and from reference 2, an esti- 

mate of the effect of the gap at M = 1.00 was made. The gap reduced 
the parameters by about 30 percent (about 15 percent lower than the 
model values). However, since the flow is viscous and, therefore, the 
effective gap is smaller than the measured gap, the tail parameters 
obtained from of reference 2 and corrected for an effective gap 

should, in general, be close to the values obtained from the model. 
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Damping in Pitch 

The variations of the amplitude ratio (the ratio of the amplitude 
from the trim line of successive oscillations to the initial oscilla- 
tion) with time, for both tail settings, are shown in figure 19. The 
time to damp to one-half amplitude, represented by the time required 
for the faired line to cross 0-50 amplitude ratio (fig. 19 ), as a 
function of Mach number, and the pitch damping-moment factors are shown 
in figure 20. As with the lift-curve slope and period data, nonlinearity 
with the trim lift coefficient was also present in the time to damp to 
one-half amplitude over all but the higher Mach numbers . Where non- 
linearity occurred, the model exhibited irregular damping character- 
istics with the higher- lift tail setting having the greater damping. 
Since, for the configuration used in this investigation, damping is 
primarily due to the tail and would be expected to decrease as downwash 
over the tail decreases, and, since, in the discussion of lift, downwash 
was believed to decrease over the tail with increasing angle of attack, 
the model damping at the higher-lift tail setting would be expected to 
be less than the damping at the low- lift tail setting. The reasons for 
the altered damping characteristics are not known at this time; however, 
since the model exhibited both longitudinal and lateral oscillations, 
the possibility exists of coupling between these oscillations. The 
effect of coupling is believed to tend to reduce the model damping. 


Directional Static Stability 

Although the primary purpose of the model flown in this investiga- 
tion was longitudinal stability, lateral oscillations were present and, 
therefore, were investigated. These oscillations are believed to be 
caused by some asymmetry in the horizontal- tail stops that could be 
present and cause lateral disturbances. 

The variation of the measured periods of the lateral oscillations 
as a function of Mach number and the static-directional-stability coeffi- 
cient (derived as in ref. J ) as determined from these periods are pre- 
sented in figures 21 and 22, respectively. Throughout the Mach number 
range the measured periods of the lateral oscillation are of the same 
order of magnitude as the measured periods of the longitudinal oscilla- 
tions and appear to be unaffected by lift coefficient over the lift 
range covered. Although some coupling between the longitudinal and 
lateral oscillations is probably present, the effect of this coupling 
is believed to be small on all the aerodynamic parameters with the 
exception of the damping in pitch. 

The static-directional-stability coefficient exhibited a decrease 
of about 30 percent from the high subsonic to the highest supersonic 
Mach numbers. Flexibility becomes increasingly important with increasing 
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Mach number, and the decrease in the static-directional-stability coef- 
ficient is believed to be due in part to flexibility of the vertical 
tail. Reference to the table of natural frequencies indicates that the 
vertical tail is rather flexible. At the Mach numbers covered, however, 
the model exhibited stable static directional characteristics. 

An estimation of the approximate maximum amplitudes of the angle 
of sideslip over the Mach number range covered has been made. At the 
high subsonic Mach numbers, maximum amplitudes of ±2° and ±3° occurred 
for the higher- and low- lift tail settings, respectively; these ampli- 
tudes decreased to about half at the transonic and supersonic Mach 
numbers . 


CONCLUSIONS 


A flight investigation of the aerodynamic characteristics at 
transonic and supersonic speeds of a rocket-propelled airplane con- 
figuration having a 52 . 5 ° delta wing and a low, swept horizontal tail 
indicated the following conclusions: 

1. At transonic speeds and with increasing Mach number, the model 
exhibited a trim change of approximately 1 ° nose-up in the low lift 
range and a 2 ° nose-down trim change while trimmed at the higher lift 
condition. 

2. The lift-curve slopes were nonlinear in the transonic and low 
supersonic Mach number range and increased with lift coefficient over 
the lift range covered . 

3. The minimum drag coefficient increased from 0.0l6 to 0.035 
through the transonic speed range. The minimum drag coefficient was 
close to the drag coefficient at zero lift. 

4. Over the Mach number range where it could be determined, the 
model exhibited poor leading-edge suction. 

5. The aerodynamic center moved rearward with increasing Mach num- 
ber up to a Mach number of about 1 . 5 , with greater stability at the 
higher lifts . 

6 . The model damping characteristics were irregular and altered 
for the two tail settings over the Mach number range covered with the 
higher-lift tail setting having the greater damping. 
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7 • The model exhibited stable static directional characteristics 
over the Mach number and lift ranges covered. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 7, 195I4.. 
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Wing 

Aspect ratio 
Area (total) 
Area (exp) 
Dihedral 
M.A.C. 

Airfoil section 


3.08 

2.976 sq ft 
I 940sq ft 
0.0 deg 
1.31 ft 
NACA 
65A 003 


Horizontal tail 
Aspect ratio 4.00 

Area (total) 0.905 sq ft 

Dihedral -20.0 deg 
M.A.C. 0.504ft 

Airfoil section NACA 

65A 006 


Vertical tail 

Area (total) 1 .37 sq ft 
Airfoil section NACA 

65A 003 


Figure 1.- General arrangement of model. All dimensions are in inches. 
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Model 



Equivalent body 



Area distribution 

Figure 2.- Nondimensional equivalent body and area distribution of model. 
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(a) Three-quarter front view. 


l-79378.1 


Figure 3«- Photographs of models. 
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(b) Model on launcher. 

Figure J.- Concluded. 
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Figure 4.- Variation of dynamic pressure with Mach number. 



Figure 5 .- Variation of Reynolds number with Mach number. 
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3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 

Time, sec 


Figure 6.- Time history during first portion of flight of some quantities 

obtained in present investigation. 
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Trim angle of attack Trim lift coefficient 
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(a) Trim lift coefficient. 
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(b) Trim angle of attack. 

Figure Longitudinal trim characteristics as a function of Mach number. 
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(a) Transonic Mach numbers. 



M » 1.23 1.29 1.34 1.41 1.47 1.55 1.65 1.76 1.33 

Angle of attack, deg 

(b) Supersonic Mach numbers. 

Figure 8.- Variation of lift coefficient with angle of attack at transonic 

and supersonic Mach numbers. 
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Figure 9*- Variation of lift- curve slope with Mach number. 
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Figure 10.- Variation of drag with lift. 
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Mach number 

(a) Maximum lift-drag ratios. 


Mach number 

(b) Lift coefficients at which maximum lift-drag ratios occur. 

Figure 11.- Variation of maximum lift-drag ratios and lift coefficients 
at which maximum lift-drag ratios occur as a function of Mach number. 
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(a) Minimum drag coefficient. 
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(b) Lift coefficient at minimum drag coefficient. 


Figure 12.- Variation of minimum drag coefficient and lift coefficient 
at minimum drag coefficient as a function of Mach number. 


CONFIDENTIAL NACA RM L54A20 


CONFIDENTIAL 


•4 















































on 

5^- 
5^ i 

5.2° 
D . 1 0 


N 

> 



- 



-B — 


: n 

— 















c 

* ' — 








0 

i — i 



















T 














vrrrqr n 


























— 

— 

1 

ttA 














































°*'Z - 8 *9 1.0 l.l 1.2 1.3 1.4 1.5 1.6 1.7 1.8 


Mach number 


Figure 13 .- Effect of lift on drag as a function of Mach number. 



Figure 14 .- Variation of measured period with Mach number. 
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Figure 15 Variation of static— stability parameter with. Mach number. 
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Figure 1 6 .- Variation of aerodynamic center with Mach number. 
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Figure 17 •- Variation of total pitching-moment coefficient with lift 

coefficient. 
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(a) Ability to produce lift. 



(b) Effectiveness in producing moment. 


Figure 18.- Variation of ability of horizontal tail to produce lift and 
effectiveness in producing moment as a function of Mach number. 
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Figure 19.- Variation of amplitude ratio with time. 
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(b) Pitch damping-moment factor. 


Figure 20.- Variation of damping characteristics of short period oscilla- 
tions with Mach number. 
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Figure 21.- Variation of measured period of lateral oscillation with 

Mach number. 



Figure 22.- Variation of static-directional-stability coefficient with 

Mach number. 
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